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Geometrical theory of whispering gallery modes 
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Abstract — Using quasiclassical approach rather precise analyt- 
ical approximations for the eigenfrequencies of whispering gallery 
modes in convex axisymmetric bodies may be found. We use the 
eikonal method to analyze the limits of precision of quasiclassical 
approximation using as a practical example spheroidal dielectric 
cavity. The series obtained for the calculation of eigenfrequencies 
is compared with the known series for dielectric sphere and with 
numerical calculations. We show how geometrical interpretation 
allows expansion of the method on arbitrary shaped axisymmetric 
bodies. 



I. Introduction 

Submillimeter size optical microspheres made of fused 
silica with whispering gallery modes (WGM)[1] can have 
extremely high quality-factor, up to 10 10 that makes them 
promising devices for applications in optoelectronics and 
experimental physics. Historically Richtmyer[2] was the first 
to suggest that whispering gallery modes in axisymmetric 
dielectric body should have very high quality-factor. He ex- 
amined the cases of sphere and torus. However only recent 
breakthroughs in technology in several labs allowed producing 
not only spherical and not only fused silica but spheroidal, 
toroidal [3], [4] or even arbitrary form axisymmetrical optical 
microcavities from crystalline materials preserving or even 
increasing high quality factor [5]. Especially interesting are 
devices manufactured of nonlinear optical crystals. Microres- 
onators of this type can be used as high-finesse cavities 
for laser stabilization, as frequency discriminators and high- 
sensitive displacement sensors, as sensors of ambient medium 
and in optoelectronical RF high-stable oscillators. (See for 
example materials of the special LEOS workshop on WGM 
microresonators [6]). 

The theory of WGMs in microspheres is well established 
and allows precise calculation of eigenmodes, radiative losses 
and field distribution both analytically and numerically. Unfor- 
tunately, the situation changes drastically even in the case of 
simplest axisymmetric geometry, different from ideal sphere 
or cylinder. No closed analytical solution can be found in this 
case. Direct numerical methods like finite elements method are 
also inefficient when the size of a cavity is several orders larger 
than the wavelength. The theory of quasiclassical methods 
of eigenfrequencies approximation starting from pioneering 
paper [7] have made a great progress lately [8]. For the 
practical evaluation of precision that these methods can in 
principal provide, we chose a practical problem of calculation 
of eigenfrequencies in dielectric spheroid and found a series 
over angular mode number I. This choice of geometry is 
convenient due to several reasons: 1) other shapes, for example 
toroids [4] may be approximated by equivalent spheroids; 2) 
the eikonal equation as well as scalar Helmholtz equation (but 
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not the vector one!) is separable in spheroidal coordinates 
that gives additional flexibility in understanding quasiclassical 
methods and comparing them with other approximations; 3) 
in the limit of zero eccentricity spheroid turns to sphere for 
which exact solution and series over I up to Z~ 8 / 3 is known 
[9]. 

The Helmholtz vector equation is unseparable [10] in 
spheroidal coordinates and no vector harmonics tangential 
to the surface of spheroid can be build. That is why there 
are no pure TE or TM modes in spheroids but only hybrid 
ones. Different methods of separation of variables (SVM) 
using series expansions with either spheroidal or spherical 
functions have been proposed [11], [12], [13]. Unfortunately 
they lead to extremely bulky infinite sets of equations which 
can be solved numerically only in simplest cases and the 
convergence is not proved. Exact characteristic equation for 
the eigenfrequencies in dielectric spheroid was suggested[14] 
without provement that if real could significantly ease the task 
of finding eigenfrequencies. However, we can not confirm this 
claim as this characteristic equation contradicts limiting cases 
with the known solutions i.e. ideal sphere and axisymmetrical 
oscillations in a spheroid with perfectly conducting walls [15]. 

Nevertheless, in case of whispering gallery modes adjacent 
to equatorial plane the energy is mostly concentrated in 
tangential or normal to the surface electric components that 
can be treated as quasi-TE or quasi-TM modes and analyzed 
with good approximation using scalar wave equations. 

II. Spheroidal coordinate system 

There are several equivalent ways to introduce prolate and 
oblate spheroidal coordinates and corresponding eigenfunc- 
tions [16], [17], [18]. The following widely used system of 
coordinates allows to analyze prolate and oblate geometries 
simultaneously: 

x = \[{e - s)(l - v 2 )} 1/2 cos(^) 
d 



(i) 



where we have introduced a sign variable s whish is equal 
to 1 for the prolate geometry with £ G [l,oo) determining 
spheroids and r\ S [—1,1] describing two-sheeted hyperboloids 
of revolution (Fig.l, right). Consequently, s = — 1 gives 
oblate spheroids for £ € [0, oo) and one-sheeted hyperboloids 
of revolution (Fig.l, right), d/2 is the semidistance between 
focal points. We are interested in the modes inside a spheroid 
adjacent to its surface in the equatorial plane. It is convenient 
to designate a semiaxis in this plane as a and in the z- 
axis of rotational symmetry of the body as b. In this case 
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d 2 /A = s(b 2 



and eccentricity e = yfl — (a/b) 2s . The 



scale factors for this system are the following: 



h„ = 
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The scalar Helmholtz differential equation is separable 



A$ + A: 2 $ = 0. 



(3) 
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$ = 0, (4) 



where c = fcd/2. The solution is $ = R m i(c, £,)S m i(c, r^e 1 ™^ 
where radial and angular functions are determined by the 
following equations: 
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i? = 0, (5) 



S 1 = 0. (6) 



Here X m i is the separation constant of the equations which 
should be independently determined and it is a function on m, I 
and c. With substitution £ = 2r/rf the first equation transforms 
to the equation for the spherical Bessel function ji(kr) if 
d/2 — ► in which case the second equation immediately 
turns to the equation for the associated Legendre polynomials 
Pmiv) with X = 1(1 + 1)- That is why spheroidal functions 
are frequently analyzed as decomposition over these spherical 
functions. 

The calculation of spheroidal functions and of A m ; is not 
a trivial task [19], [20]. The approximation of spheroidal 
functions and their zeroes may seem more straightforward for 
the calculation of eigenfrequencies of spheroids, however we 
found that another approach that we develop below gives better 
results and may be easily generalized to other geometries. 

III. ElKONAL APPROXIMATION IN SPHEROID 

The eikonal approximation is a powerful method for solving 
optical problems in inhomogeneous media where the scale 
of the variations is much larger than the wavelength. It was 
shown by Keller and Rubinow [7] that it can also be applied to 
eigenfrequency problems and that it has very clear quasiclas- 
sical ray interpretation. It is important that this quasiclassical 
ray interpretation requiring simple calculation of the ray paths 
along the geodesic surfaces and application of phase equality 
(quantum) conditions gives precisely the same equations as the 
eikonal equations. Eikonal equations allow, however, to obtain 
more easily not only eigenfrequencies but field distribution 
also. 



In the eikonal approximation the solution of the Helmholtz 
scalar equation is found as superposition of straight rays: 



u(r) = A(r)e lkoS{r \ 



(7) 



The first order approximation for the phase function S called 
eikonal is determined by the following equation. 



(8) 



where e is optical susceptibility. For our problem of searching 
for eigenfrequencies e does not depend on coordinates, e = n 2 
inside the cavity and e = 1 - outside. Though the eikonal can 
be found as complex rays in the external area and stitched on 
the boundary as well as ray method of Keller and Rubinow 
[7], [8], [21] can be extended for whispering gallery modes 
in dielectrical bodies in a more rigorous way [22]. To do so 
we must account for an additional phase shift on the dielectric 
boundary. Fresnel amplitude coefficient of reflection [23]: 
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where \ — 1 f° r quasi-TE modes and \ = 1/n 2 for quasi-TM 
modes, give the following approximations for the phase shift 
for grazing angles: 
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However direct use of this phase shift in the equations for 
internal rays as suggestes in [22] leads to incorrect results. 
The reason is a well known Goos-Hanchen effect - the shift 
of the reflected beam along the surface. The beams behave 
as if they are reflected from the surface hold away from the 
real boundary at a r = 2fc cos e ■ ^ at ^ s wn y we ma y substitute 
the problem for a dielectric body with the problem for an 
equivalent body enlarged on o r with the totally reflecting 
boundaries. The parameters of equivalent spheroid are marked 
below with overbars. 

The eikonal equation separates in spheroidal coordinates if 
we choose S - S&) + S 2 (v) + S 3 (4>): 
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After immediate separation of ( ^ 



H we have: 
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Fig. 1. Graphical representation of prolate and oblate spheroidal coordinate systems (£, rj, <f>). 



Introducing another separation constant v we finally obtain 
solutions: 

,2 \ V 2 



S/J, 



dS 2 (r)) 



v l n 2 d 2 sri 2 



drj \1-T] 2 4(1 — 77 2 ) (l-?y 2 ) 
which after some manipulations transform to: 
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where 



(1 + sa) - + sa) 2 - 4scm7o 
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where a 



4v 2 



V 2 o 



1 — fj, 2 /v 2 . It is now the time to turn to 
the quasiclassical ray interpretation [7], [8]. The equation for 
the eikonal describes the rays that can spread inside spheroid 
along the straight line. These are the rays that freely go inside 
spheroid than touch the surface and reflect. For the whispering 
gallery modes the angle of reflection is close to ir/2. The 
closest to the center points of these rays form the caustic 
surface which is the ellipsoid determined by a parameter £ r . 
The rays are the tangents to this internal ellipsoid and follow 
along the geodesic lines on it. In case of ideal sphere all the 
rays of the same family lie in the same plane. However, even 
a slightest eccentricity removes this degeneracy and inclined 
closed circular modes which should be more accurately called 
quasimodes [24] are turned into open-ended helices winding 



up on caustic spheroid precessing [25], and filling up the whole 
region as in a clew. The upper and lower points of these 
trajectories determine other caustic surface with a parameter 
r) r determining two-sheeted hyperbolid for prolate or one- 
sheeted hyperboloid for oblate spheroid. The value of r/ r has 
very simple mechanical interpretation. The rays in the eikonal 
approximation are equivalent to the trajectories of a point-like 
billiard ball inside the cavity. As axisymmetrical surface can 
not change the angular momentum related to the z axis, it 
should be conserved as well as the kinetic energy (velocity). 
That is why r\ T is simply equal to the sine of the angle between 
the equatorial plane and the trajectory crossing the equator and 
at the same time it determines the maximum elongation of the 
trajectory from the equator plane. If all the rays touch the 
caustic or boundary surface with phases that form stationary 
distribution (that means that the phase difference along any 
closed curve on them is equal to integer times 2tt), then the 
eigenfunction and hence eigenfrequency is found. 

To find the circular integrals of phases kS dl4> we should 
take into account the properties of phase evolutions on caustic 
and reflective boundary. Every touching of caustic adds ir/2 
(see for example [8]) and reflection adds ir. Thus for Si we 
have one caustic shift of ir/2 at £ r and one reflection from 
the equivalent boundary surface £ s (at the distance a from the 
real surface), for S 2 - two times tt/2 due to caustic shifts at 
±r/ r , and we should add nothing for S3. 



fcASi = 2k 



</r 



dSi 



= 2tt(<7 - 1/4) 



kAS 2 = 2k J ^j-drj = 2tt(p + 1/2) 
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fcA& = k 



dS* 



= 2w\m\ 



(16) 



where q — 1,2,3... - is the order of the mode, showing 
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the number of the zero of the radial function on the surface, 
and p — I — |m| = 0, 1,2.... These conditions plus integrals 
(114-1 completely coincide with those obtained by Bykov [26], 
[27], [28] if we transform ellipsoidal to spheroidal coordinates, 
and have clear geometrical interpretation. The integral for Si 
corresponds to the difference in lengths of the two geodesic 
curves on r\ T between two points P\ = i] r , tfii ) and 
P*z = (6r, Vr, ^2 )■ The first one goes from the caustic circle 
of intersection between £ r and r/ r along r\ r to the boundary 
surface £ s , reflects from it, and returns back to the same circle. 
The second is simply the arc of the circle between P\ and 
The integral for 52 corresponds to the length of a geodesic 
line going from P\ along £ r , lowering to — rj r and returning 
to r\ r at P2 minus the length of the arc of the circle between 
Pi and P2 . The third integral is simply the length of the circle 
of intersection of £ r and r\ r . 

These are elliptic integrals. For the whispering gallery 
modes when r\ r <C 1 and £0 — £r "C £r> S2 may be expanded 
into series over r\ r and ( and integrated with the substitutions 
of 77 = r] r s'mip, £ = (£ 2 — £ 2 )/£ 2 an d Finally, expressing 
spheroidal coordinates £ r and expressing £0 through parame- 
ters of spheroid, we have: 
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4 + 64^ Vr 
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S 3 =n<j>, 



1536$ 



Now we should solve the following system of equations: 



A;AS' 1 =2fcS , 1 (Co) 
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3a 3 VTTCo 
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Co 372 ( 1 
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26 2 



--2ir(p + 1/2) 
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\/T+Co" 



=2ttM 



(18) 



Using the method of sequential iterations, starting for example 



from nk^a = I, Q 



(0)' 



0, r?.. 



(0) 



this system may be 



-x-FEM 

-O— Equation (19) 

-A— References [5,26-18] 
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Fig. 2. Comparison of the precision of calculation of eigenfrequencies in 
spheroid. 
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nka=nk(a — a r ) = I — (3 q 



I \ 1/3 2p(a -b) + a 



2b 



Xn 



(17) 



Vn 1 



30? 
20 V 2 



-1/3 



Pq (2p(a 3 -b 3 )+a 3 f 2n X (2 X 2 -3n 2 )\ fiy 2/3 
'12 



b 3 



{n 2 - l) 3 / 2 



-o(r 



(19) 



where for the convenience of comparison we introduced j3 q = 
— [|7r(g— -|)] 2 ^ 3 . The value of cosf? needed for the calculation 
of r (II It one may estimate as cosf? = — I 2 /(nka) ~ 
/31- 1 / 3 . 

The first three terms for nka were obtained in [3], [26], [27], 
[28] from different considerations, the last three are new. 

To test this series we calculated using finite element method 
eigenfrequencies of TE modes in spheroids with different 
eccentricities with totally reflecting boundaries for I = m = 
100 (Fig. 2). Significant improvement of our series is evident. 
The divergence of the series for large eccenricities is explained 
by the fact that the approximation that we used to calculate 
the integrals d!6i but not the method itself breaks down in this 
case. Namely i] r becomes comparable to £ and should not be 
treated as a small parameter. 

If we put a — b, then all six terms in the obtained series 
coincide with that obtained in [9] from exact solution in 
sphere with two minor differences: 1) Airy function zeroes 
a q ~ (-2.3381,-4.0879,-5.5206,...) stand in Schiller's 
series instead of approximate f3 q values (a q — (3 q ~ -0.017; 
-0.0061; -0.0033, ...). The reason is that the eikonal approx- 
imation breaks down on caustic, where its more accurate 
extension with Airy functions is finite [7]. To make our 
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solution even better we may just formally use a q instead of 
j3 q \ 2) Minor difference in the last term is caused, we think, 
by misprint in [9], where in our designations instead of e?2 = 
-2 2 / 3 n 2 (-3 + 2x 2 )a g /6 should be d 2 = -2 2 / 3 (-3n 2 + 2x 2 ). 
The eikonal equation for the sphere may be solved explicitly 
and the expansion of the solution shows that quasiclassical 
approximation breaks down on a term 0(Z _1 ), and of the same 
order should be the error introduced by substitution of vector 
equations by a scalar ones. 

It is interesting to note that when a = 2b (oblate spheroid 
with the eccentricity e = V0.75, the eigenfrequency separation 
in the first order of approximation between modes with the 
same I determined by the third term becomes equal to the 
separation between modes with different I and the same / — m 
(free spectral range). The difference appears only in the term 
proportional to 0(7~ 2 / 3 ). This situation is close to the case 
that was experimentally observed in [3]. This new degeneracy 
has simple quasigeometrical interpretation - like in case of a 
sphere geodesic lines inclined to the equator plane on such 
spheroid are closed curves returning at the same point of 
the equator after the whole revolution, crossing, however, the 
equator not twice as big circles on a sphere but four times. 

IV. Arbitrary convex body of revolution 

To find eigenfrequencies of whispering gallery modes in 
arbitrary body of revolution one may use directly the results of 
the previous section by fitting the shape of the body in convex 
equatorial area by equivalent spheroid. In fact the body should 
be convex only in the vicinity of WG mode itself. For example 
a torus with a circle of radius ro with its center at a distance 
Ro from the z axis as a generatrix may be approximated 
by a spheroid with a = Rq + ro and b = J (Rq + ro)ro- 
Nevertheless, more rigorous approach may be developed. 

The first step is to find families of caustic surfaces. This is 
not a trivial task in general but it is equivalent to finding caustic 
curves for the plane curve forming the body of revolution 
which is in fact the so-called biinvolute curve (the difference 
in length between a sum of two tangent lines from a point on 
a curve to a biinvolute curve and an arc between these lines 
is constant). Unfortunately we can not give now the formal 
provement of this statement but it looks true. Another family 
form curves orthogonal to the first family. For example in case 
of torus these families are concentric circles and radii, that is 
why caustic surfaces for a torus are also concentric toruses. 
In general case the following approximation may be used to 
find the first family of biinvolute curves [8]: 



n(s) ~ k P 1 ^(s) + 0{k 2 ), 



(20) 



where n(s) is the normal distance from a point s on the curve 
to a biinvolute curve, k is a parameter of a family, and p is 
the radius of curvature of the initial curve at s. 

Let we have found a caustic surface from the first family, 
parametrized as 



z = u 

x = g(u) cos 4> 

y = g(u) sin (f> 



(21) 



A geodesic line for this surface is given by the following 
integral: 



cr 



y/T+g 3 



(22) 



du g(u)^/g 2 {u)-c 2 

where c\ is some constant, which is equal in our case to p r = 
g(zmax) - the radius of caustic circle at maximum distance 
from equatorial plane. The length of geodesic line: 




vV - Pr 

The length of geodesic line, connecting points <j>\ and $2 ■ 



L\ = 2 



VV - Pr 



du 



(24) 



The length of arc from <po — to <p c — 2 f^ r -^du is equal 



to L% = p r (j) c . 



Finally: 



Lf = 2 



9\j9 2 - Pr 



du 



(25) 



nk(Lf 



Ink 



V^ + g'Wd 2 - A 



du 



2ir(p + 1/2) 



(26) 



In analogous way for another geodesic line on a caustic 
surface from the other family: 



z = v 

x = h(v)coscj) 
y = h(v)sm<j), 



(27) 



we have 

nfc(Lj - L\) 



2nk VTTk^ER dv 



= 277(9 -1/4) 



(28) 



The third condition is: 



2imkp r = 27r|m| (29) 
With the substitution u = j^r] r , v = f£ r ?7, <?(«) = 

u/e^ ~s^r^w r , h{u) = i^—s^T^f in Gsm 

U9\ we again obtain expressions for spheroid obtained before. 
For a torus caustic surfaces are toruses and cones and are 
determined by the equations: 

g[u) = Ro + \Jr 2 -u 2 

Pr — Ro 



h(v) = Ro + v- 



(30) 



Vr 2 o - (Pr - Ro) 2 
In conclusion. We have analyzed quasiclassical method 
of calculation of eigenfrequencies in spheroidal cavities and 
found that it gives approximations correct up to the term 
proportional to l~ 2 ^ 3 . This method may be easily expanded 
on arbitrary convex bodies of revolution. 



Acknowledgment 

The work of M.L.Gorodetsky was supported in part by the 
Alexander von Humboldt foundation. 

References 

[1] V. B. Braginsky, M. L. Gorodetsky and V. S. Ilchenko, "Quality-factor 
and nonlinear properties of optical whispering-gallery modes," Phys. 
Lett. A137, pp. 393-397, 1989. 

[2] R.D. Richtmyer, "Dielectric Resonators", J. of Appl. Phys. 10, pp. 391— 
398, 1939. 

[3] V.S. Ilchenko, M.L.Gorodetsky, X.S.Yao and L.Maleki, "Microtorus: a 

high-finesse microcavity with whispering-gallery modes", Opt. Lett. 26, 

pp. 256-258, 2001. 
[4] K. Vahala, "Optical microcavities", Nature 424, pp. 839-846, 2001. 
[5] V.S. Ilchenko, A.A. Savchenkov, A.B. Matsko et al. "Nonlinear optics 

and crystalline whispering gallery mode cavities," Phys. Rev. Lett. 92, 

(043903), 2004. 

[6] 2004 Digest of the LEOS Summer Topical Meetings: Biophoton- 
ics/Optical Interconnects & VLSI PhotonicsAVGM Microcavities (IEEE 
Cat. No.04TH8728), 2004. 
[7] J.B. Keller, S.I. Rubinow, "Asymptotic solution of eigenvalue problems", 

Ann. Phys. 9, pp. 24-75, 1960. 
[8] V.M. Babic, V.S. Buldyrev, Short-wavelength diffraction theory. Asymp- 
totic methods, Springer- Verlag, Berlin Heidelberg, 1991. 
[9] S.Schiller, "Asymptotic expansion of morphological resonance frequen- 
cies in Mie scatternig", Appl. Opt. ,32, pp. 2181-2185, 1993. 

[10] R. Janaswamy,"A note on the TE7TM decomposition of electromagnetic 
fields in three dimensional homogeneous space", IEEE Trans. Antennas 
and Propagation 52, pp. 2474-2477, 2004. 

[11] S.Asano, G.Yamamoto, "Light scattering by a spheroidal particle", Appl. 
Opt. 14, pp. 29^19, 1975. 

[12] VG. Farafonov, N.V Voshchinnikov, "Optical properties of spheroidal 
particles", Astrophys. and Space Sci. 204, pp. 19-86, 1993. 

[13] A.Charalambopoulos, D.I.Fotiadis, C.V Massalas, "On the solution of 
boundary value problems using spheroidal eigenvectors", Comput. Phys. 
Comm. 139, pp. 153-171, 2001. 

[14] RC.G. de Moraes, L.G. Guimaraes, "Semiclassical theory to optical 
resonant modes of a transparent dielectric spheroidal cavity", Appl. Opt. 

41, pp. 2955-2961, 2002. 

[15] L. Li, Z. Li, M. Leong, "Closed-form eigenfrequencies in prolate 

spheroidal conducting cavity", IEEE Trans. Microwave Theory Tech. 

51, pp. 922-927, 2003. 
[16] I. V. Komarov, L. I. Ponomarev, and S. J. Slavianov, Spheroidal and 

Coulomb Spheroidal! 'unctions (in russian) (Moscow: ), Nauka, Moscow, 

1976. 

[17] L.Li, X.Kang, M.Leong, Spheroidal Wave Functions in Electromagnetic 

Theory, John Wiley & Sons, 2002. 
[18] Handbook of Mathematical Functions, ed. M.Abramowitz and 

I.E.Stegun, National Bureau of Standards, 1964. 
[19] RC.G. de Moraes, L.G. Guimaraes, "Uniform asymptotic formulae for 

the spheroidal radial function", /. of Quantitative Spectroscopy and 

Radiative Transfer 79-80, pp. 973-981, 2003. 
[20] RC.G. de Moraes, L.G. Guimaraes, "Uniform asymptotic formulae for 

the spheroidal angular function", /. of Quantitative Spectroscopy and 

Radiative Transfer 74, pp. 757-765, 2003. 
[21] V.A.Borovikov, B.E.Kinber, Geometrical theory of diffraction, IEE Elec- 
tromagnet. Waves Ser.37, London, 1994. 
[22] E.L. Silakov,"On the application of ray method for the calculation of 

complex eigenvalues" (in russian), Zapiski nauchnogo seminara LOMI 

42, pp. 228-235, 1974. 

[23] M.Born and E.Wolf, Principles of Optics, 7-th ed., Cambridge University 
Press, 1999. 

[24] VI. Arnold, "Modes and quasimodes" (in russian), Funktsionalny analiz 

i prilozheniya 6, pp. 12-20, 1972. 
[25] M.L.Gorodetsky, V.S. Ilchenko, "High-Q optical whispering-gallery mi- 

croresonators: precession approach for spherical mode analysis and 

emission patterns with prism couplers", Opt. Commun. 113, pp. 133- 

143, 1994. 

[26] V.P.Bykov, "Geometrical optics of three-dimensional oscillations in open 
resonators" (in russian), Elektronika hol'schikh moschnostei 4, pp. 66- 
92, Moscow, 1965. 

[27] L.A. Vainstein,"Barellike open resonators" (in russian), Elektronika 
hol'schikh moschnostei 3, pp. 176-215, Moscow, 1964. 

[28] L.A. Vainstein, Open Resonators and Open Waveguides, Golem, Denver, 
1969. 



